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> ■ We present seven epochs of spectroscopy on the quadruply imaged quasar 
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SDSS J1004+4112, spanning observed-frame time delays from 1 to 322 days. 



The spectra reveal differences in the emission lines between the lensed images. 
Specifically, component A showed a strong enhancement in the blue wings of 
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several high-ionization lines relative to component B, which lasted at least 28 
days (observed frame) then faded. Since the predicted time delay between A and 
B is <30 days, our time coverage suggests that the event was not intrinsic to 
the quasar. We attribute these variations to microlensing of part of the broad 
emission line region of the quasar, apparently resolving structure in the source 
plane on a scale of ~10 16 cm at z — 1.734. In addition, we observed smaller dif- 
ferences in the emission line profiles between components A and B that persisted 
throughout the time span, which may also be due to microlensing or millilensing. 
Further spectroscopic monitoring of this system holds considerable promise for 
resolving the structure of the broad emission line region in quasars. 

Subject headings: gravitational lensing — quasars: general — quasars: emission 
lines — quasars: individual (SDSS J100434.91+41 1242.8) 

1. Introduction 

Microlensing in the images of a multiply-imaged quasar was first reported by Irwin et al. 
(1989) for the quadruple lens Q2237+0305 (Huchra et al. 1985). Most quasar microlensing 
studies have been based on broad-band photometric monitoring (e.g., Wozniak et al. 2000; 
Schechter et al. 2003; Colley & Schild 2003), which is sensitive primarily to variations in 
the continuum. Microlensing of the continuum is expected since the optical/UV continuum 
emission is thought to originate in a region that is comparable in size to the Einstein radius 
of a typical star in a typical lens galaxy. 

Microlensing of the broad emission line region (BELR) is also possible, if the BELR has 
structure on scales comparable to the Einstein radius of a star (Nemiroff 1988; Schneider 
& Wambsganss 1990). The possibility of BELR microlensing seemed rather remote until 
recent reverberation mapping work revised the estimate of the BELR size downward from 
~10 18 cm to ~10 16 cm (Wandel, Peterson, & Malkan 1999; Kaspi et al. 2000). Inspired by 
these numbers, Abajas et al. (2002) and Lewis & Ibata (2003) revived the idea of looking for 
microlensing of the BELR and computed possible line profile variations for various BELR 
models. 

Possible examples of microlensing of a quasar emission line have been presented by Fil- 
ippenko (1989) for Q2237+0305, and by Chartas et al. (2003) for H1413+117. In particular, 
Chartas et al. (2003) detected a strong, redshifted Fe Ka emission line in the X-ray spectrum 
of only one of the components in the quadruple lens H1413+117. Although they did not 
have multiple epochs to look for variability, Chartas et al. (2003) invoked a short predicted 
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time delay between the components to argue that microlensing is the preferred explanation 
for seeing the Fe Ka line in only one component. 

In this paper we present results from spectroscopic monitoring of the recently discovered 
quadruple lens SDSS J1004+4112 (Inada et al. 2003; Oguri et al. 2003). We observe 
variability in the broad emission line profiles of one of the lensed images that provides 
strong evidence for microlensing of the BELR, suggesting that the theoretical predictions for 
microlensing were correct and confirming that the BELR has structure on the scale of the 
Einstein radius of a star. 



2. The Data and Spectral Analysis 

Component B of SDSS J1004+4112 was first identified as a quasar by Cao, Wei, & 
Hu (1999) and was also targeted as a quasar candidate (Richards et al. 2002) as part of 
the Sloan Digital Sky Survey (SDSS; York et al. 2000); the SDSS spectrum was taken 
on 2003 February 3 and the object was identified as a z = 1.734 quasar. Component A 
was also identified as a possible quasar candidate based on its colors in the SDSS imaging 
data (Fukugita et al. 1996; Gunn et al. 1998); it was first confirmed as a quasar (also 
at z — 1.734) from observations taken on 2003 May 3 using the ARC 3.5m telescope at 
Apache Point Observatory. Higher quality spectra of each of the four components and the 
lens galaxy were taken on 2003 May 31 with the LRIS (Oke et al. 1995) spectrograph on 
the Keck I telescope at the W. M. Keck Observatory. These spectra are further described 
in Inada et al. (2003), Oguri et al. (2003), and Table 1. Table 1 gives the UT date of the 
observations, the telescope used, the components observed, the exposure time, the spectral 
dispersion, and the wavelength range covered. 

Inada et al. (2003) and Oguri et al. (2003) noted that the Keck/LRIS spectra showed 
that all four components had the same redshift and similar spectra, but there are some 
obvious differences — with component A showing the largest differences relative to the other 
components (see Fig. 1). Specifically, component A has a much stronger blue emission line 
wing in the high-ionization lines (Si IV/O IV], C IV, and He II). Component B appears 
to have a slightly enhanced red wing as compared with the other three components. In 
Figure 1, we have sought to emphasize the differences in the C IV emission line profile by 
subtracting a continuum fit between A 1450 A and A 1690 A, and then normalizing to the 
peak of C IV emission. Other choices in presentation could enhance or reduce the appearance 
of similarities and differences between the components. 

Based on the puzzling differences in the broad emission line profiles of the four compo- 
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nents of the Keck/LRIS spectra, in 2003 November and 2003 December we obtained addi- 
tional spectra of images A and B to monitor the differences; see Table 1. We used the DIS 111 
spectrograph on the ARC 3.5m telescope at Apache Point Observatory, using the same in- 
strument setup and data reduction for each epoch. The dispersion was 2.4 A pixel and the 
spectra cover the range AA3890-9350 A. Flux calibration was performed with respect to a 
hot standard star: either Feige 34 or G 191-B2B. Wavelength calibration was performed with 
respect to a helium-neon-argon comparison lamp. Both components (which are separated by 
3'/73) were observed at the same time in order to minimize any differences between them, at 
the expense of loss of spectrophotometric accuracy (slit not positioned at the parallactic an- 
gle). At each epoch we took two spectra of either 40 or 45 minutes each. These spectra were 
extracted separately and then combined to yield 2003 November /December epoch spectra 
shown in Figure 2. The 2D spectra are only moderately resolved, so we deblended the 2D 
spectra by fitting a double Gaussian profile before extracting the spectra. We estimate that 
contamination of the faint component by the brighter component is less than 3%. 

Figure 2 shows all seven epochs of spectroscopic data for components A and B in the C IV 
emission line region. The 2003 May 3 epoch observation of component A at APO confirms 
the reality of the excess emission in the blue wing of component A in the Keck spectrum. 
The excess emission in He II spans a velocity range ~2500-8500 km s _1 (blueshifted) . Both 
Si IV and C IV have excesses over a slightly higher velocity range. The excess in He II is 
much stronger in terms of equivalent width than that in C IV. 

3. Discussion 

3.1. Time-dependent Blue Wing Differences: Microlensing of the BELR 

3.1.1. Discounting Alternatives and Objections 

To argue that the excess in the blue wings of the high-ionization emission lines of com- 
ponent A is caused by microlensing of the BELR, we must rule out alternative explanations. 
The first question is whether the flux enhancements are an artifact of the data acquisition 
or reduction procedures. We view this explanation as highly unlikely because the enhance- 
ment appeared at two different epochs, in spectra taken by different observers with different 
telescopes, and reduced independently. 

The second question is whether the spectral variations could be intrinsic to the quasar, 
rather than induced by lensing. If so, then the same variations should be seen in the other 
lensed images, offset in time by the lens time delays. While the time delays are still uncertain, 
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in nearly all of the lens models presented by Oguri et al. (2003) the delay between components 
A and B is predicted to be <30 days. (The Oguri et al. models do not form an exhaustive 
set, but the prediction of a short time delay between A and B is generic.) In ~10% of the 
models component A leads component B. In this case seeing the flux enhancement in the 
2003 May 3 and 2003 May 31 spectra of component A but not in the 2003 May 31 spectrum 
of component B means that intrinsic variability cannot explain the data. In the other 
~90% of the models component B leads component A. In this case, the intrinsic variability 
hypothesis would imply that the variations must have appeared in component B sometime 
after 2003 February 3, been present during 2003 April (in advance of the variations observed 
in component A in 2003 May), and then disappeared before 2003 May 31. In other words, 
the event must have lasted <117 days in the observed frame, or <43 days in the rest frame 
- and we were fortunate to catch the event in component A just before it disappeared. 

Since the time separation of our spectral coverage does not completely rule out intrinsic 
variability, this issue deserves further discussion. Specifically, we must address the likelihood 
of intrinsic variability of the BELR on timescales of less than 43 days in the rest frame. 
Reverberation timescales can be less than 43 days since the radius of the BELR is likely to 
be on the order of or less than this size. However, in that case the entire emission line profile 
will vary with respect to any significant change in the continuum. Since we do not observe 
a significant change in the continuum level (see below) and since the enhancement of the 
emission lines is only in the blue wing and not over the entire profile, a reverberation effect 
is unlikely. Thus we are left with the possibility of a dynamical change — such as is seen in 
the so-called double-peaked emission line quasars (e.g., Eracleous & Halpern 2003; Strateva 
et al. 2003) — which could better explain the blue-wing-only nature of the enhancement in 
the high-ionization emission lines. Indeed the dynamical timescale (~ 6 months; Eracleous 
2003) can be in the range needed to explain our observations. However, double-peaked 
emission is typically only seen in the Balmer lines (and sometimes Mg II), and is absent or 
weak at best in high-ionization lines like C IV and He II. Furthermore, it would be very 
unusual for a double-peaked emission line object to show a blue-wing enhancement in both 
our 2003 May 3 and 2003 May 31 epochs, but no strong blue- or red-wing enhancement of 
our four 2003 November/December epochs. Thus, although an explanation of the data in 
terms of intrinsic variability is not rigorously excluded, such a model would have to be rather 
improbably contrived. 

One possible objection to the microlensing hypothesis is if the BELR was microlensed, 
why not the continuum as well? It is difficult to put constraints on any enhancement of 
the continuum of component A relative to B during the time spanned by our observations, 
because not all of the spectra were taken at the parallactic angle. Still, we can estimate 
that component A was no more than ~20% brighter than B in the continuum in 2003 May 
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(as compared to 2003 November/December). This is not obviously inconsistent with the 
microlensing hypothesis, however. It is easy to imagine configurations in which part of 
the BELR is close enough to the caustic in the source plane of a star in the lens to be 
micro lensed, while the continuum source was far enough to feel little effect. In fact, this 
picture is consistent with our hypothesis that only part of the BELR was being microlensed 
in 2003 May (see § 3.1.3). 

Another possible complication with the microlensing hypothesis is that cluster galaxies 
tend to be stripped of their halos through mergers and interactions. Microlensing would 
therefore require either a cluster member very close to the line of sight to component A or 
a population of intracluster stars or massive compact halo objects (MACHOs; e.g., Totani 
2003; Baltz et al. 2003). The presence of intracluster MACHOs might not be surprising, 
because tidal forces would naturally strip MACHOs from cluster galaxies along with the 
galaxy halos. On the other hand, there is evidence for a galaxy superimposed on component 
A that could host the microlensing object, as shown in Figure 3. Proper PSF subtraction to 
confirm this hypothesis is not possible since component A is saturated in this image, but our 
best efforts do reveal residual flux at the center of the circle in Figure 3 with an estimated 
magnitude of i ~ 24.5. 



We can also reverse the argument and ask if our microlensing hypothesis makes sense 
given our current understanding of the structure of quasars and the details of microlensing. 
Without invoking detailed microlensing and BELR models (which would be premature since 
we have observed only one microlensing event), we still find that qualitative and quantitative 
arguments reveal good consistency between the microlensing hypothesis and the data. 



The scale for microlensing is given by the Einstein radius of a star (Schneider, Ehlers, 
& Falco 1992), 



where -REin,L and -REin,s are the Einstein radius projected into the lens and source planes, 
respectively, and D h D s , and D is are the angular diameter distances to the lens, to the 
source, and from the lens to the source, respectively. Microlensing is said to occur at low 
optical depth if the mean separation between stars is d ^> -REin,L and stars contribute a 



3.1.2. Examining the Microlensing Hypothesis 
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small fraction of the surface mass density, or at high optical depth if d ~ -REin,L and stars 
contribute a substantial fraction of the surface mass density. The microlensing probability 
can be high even if the optical depth is low (Schechter & Wambsganss 2002). Regardless of 
whether the optical depth is high or low, the Einstein radius sets the characteristic scale, 
and the BELR can be microlensed if it has structure on scales <-REin,s (Abajas et al. 2002; 
Lewis & Ibata 2003). For SDSS J1004+4112, the lens redshift is z x = 0.68 and the source 
redshift is z s = 1.734, so the Einstein radius of a 0.1 M star associated with the lensing 
cluster is 

OlMoJ k ™ /2 Cm ~ 5 ' 3 (olMo ) ^o /2 It-days (3) 

in a cosmology with Qm = 0.3, Q\ = 0.7, and Hq = 70 /i7okms _1 Mpc -1 . Thus, the typical 
stellar Einstein radius is indeed comparable to the currently favored size of the BELR, and 
microlensing is not unexpected. 

The variability is caused by relative motion between the caustic network and the source. 
If the motion is dominated by the proper motion of the lens galaxy (with transverse velocity 
v±), then the effective transverse velocity projected into the source plane and expressed in 
distance per unit observed-frame time is 

= ~r+Ti ~Di ' w 

If the microlensing is dominated by caustic crossings, then the characteristic event duration 
is the time for the caustic to sweep across the source, t src ~ 2R src /v e s. If we estimate v ± ~ a 
where a is the velocity dispersion of the lens, then a ~ 700 km s" 1 for the lensing cluster in 
SDSS J1004+4112 (Oguri et al. 2003) yields 

Ihn \ ( V± 

10 16 cm I \ 700 km s" 



*- ~ I T7g^ I I ,„„ ;i -, I 3- (5) 



The fact that we see variability in the broad emission lines on a time scale of ~6 months 
suggests either that we happened to catch the end of a long event, or that microlensing is 
affecting only part of the BELR with a characteristic size smaller than 10 16 cm (or both). 

A second interesting time scale is the typical time between microlensing events. Here 
the key unit is the time to cross an Einstein radius, 

o / oaf M V /2 f V ± Y\-i /2 fR . 

^Ein ~ -#Ein s Vett ~ 8.6 — — T h 7n ' yr 6 

\0.1M Q ) V700 kms~V 70 KJ 

If even a few percent of the mass is in stars, microlensing is associated with a caustic 
network rather than a single star, and so naive estimates of the time between events are 
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difficult. To obtain a better estimate, we have used the standard ray shooting technique to 
compute microlensing magnification maps and generate sample light curves (Kayser, Refsdal, 
& Stabell 1986; Wambsganss, Paczynski, & Katz 1990). We find that starting from a non- 
microlensed position, the average wait time until the magnification changes by >30% is 
~ (0.2-0.8) x t-E m - The range represents uncertainties in the source size and in the number 
density of microlenses. Although there are many uncertainties, it seems reasonable to expect 
that microlensing events could be observed in SDSS J1004+4112 on an approximately yearly 
basis. Incidentally, if t src > t Ein then microlensing events will blur together and the light 
curves will show continuous smooth variations (see, e.g., Kochanek 2003). The apparent lack 
of microlensing in the 2003 November/December epochs therefore adds further support to 
the hypothesis that the region being microlensed is smaller than ~10 16 cm. 

The microlensing hypothesis nicely explains one of the more interesting observational 
results, namely that different emission lines show different amounts of variability. From 
reverberation mapping results, we know that the BELR is stratified by ionization and that 
higher ionization lines are found closer to the center (Peterson & Wandel 1999). As a 
result, the highest ionization regions have the smallest effective sizes and should be the most 
sensitive to microlensing. Our finding that the excess in the He II line is stronger than 
the excess in the C IV line (in terms of equivalent width) is consistent with their relative 
reverberation mapping sizes (Peterson & Wandel 1999). Furthermore, there is a suggestion 
of a weaker excess in the blue wings of the C III] /Si III]/A1 III/FeIIIUV34 complex, Mg II, 
and Fe III UV48 (a triplet at AAA 2062.211, 2068.904, 2079.652); see Figure 4, and also the 
ratio spectra in Figure 3 of Oguri et al. (2003). All the above lines are lower ionization lines 
than C IV and He II, so their weaker excesses are also consistent with being more weakly 
microlensed as a result of their emitting regions being larger with respect to the projected 
Eistein radius of the lens. 

An obvious question is why should microlensing of the BELR be seen in this lens system, 
but not in others? One possibility has to do with the quasar's intrinsic luminosity. The 
observed absolute i magnitude of component A is Mi = —26.9, but this component is 
amplified by a factor of ~20 or more, 1 so the intrinsic absolute i magnitude is Mj > —23.8. 
Many other lensed quasars have similar observed magnitudes but smaller amplifications, 
and hence higher luminosities. Since the size of the BELR scales as the 0.5-0.7 power of 
the luminosity (e.g., Kaspi et al. 2000), SDSS J 1004+41 12 being relatively under- luminous 



1 An isothermal ellipsoid plus shear lens model, which is simple but fits the data well, implies an amplifi- 
cation factor of 23 for component A. More complicated models yield a broad range of amplifications where 
the median is 18 but there is a long tail to amplifications of 100 or more. See Oguri et al. (2003) for modeling 
details. 
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would make it more sensitive to microlensing than many other lensed quasars. 

The large amplification helps make microlensing unsurprising in another way. Large 
amplifications are associated with significant distortions that increase the size of the caustics 
(see Fig. 2 of Schechter & Wambsganss 2002). Thus, with a given star field the microlens- 
ing probability increases as the amplification increases. This may explain not only why 
microlensing has been detected in SDSS J1004+4112, but also why it was seen first in com- 
ponent A (the highest-amplification component). 

3.1.3. Implications for BELR Structure 

Obtaining concrete constraints on the structure of the BELR will require detailed mod- 
eling, and would greatly benefit from observations of additional microlensing events. Nev- 
ertheless, combining the data in hand with general arguments already permits some strong 
and valuable conclusions. First, unless the number of clouds is orders of magnitude smaller 
than required by other high- resolution Keck observations (Arav et al. 1998), the asymmetric 
nature of the microlensing (lensing of the blue wing only) rules out a spherically symmet- 
ric distribution of dynamically virialized, thermal line-width clouds (Nemiroff 1988; Abajas 
et al. 2002; Lewis & Ibata 2003). 

Pure radial outflow models can be ruled out because they would produce symmetric line 
profile changes with most of the variation in the line core (Nemiroff 1988). Pure radial inflow 
models cannot be fully excluded, but asymmetric microlensing would require the accretion 
disk to have a radial extent comparable to the radial extent of the BELR and to be optically 
thick (such that we do not see clouds at all velocities). 

The asymmetry therefore seems to imply that a strong rotational component is needed 
in the high-ionization region of the BELR. Such a component could come in the form of 
a pure Keplerian disk or a rotating disk- wind (e.g., Murray & Chiang 1998; Elvis 2000). 
Microlensing of the part of the BELR that is rotating toward us would then easily explain 
the features that we observe. 

The most robust statement we can make is that the observations confirm that the 
BELR has structure on the scale of the Einstein radius. Because we are sampling a region 
on the order of ~ 1.4 x 10 16 cm, we effectively have a "telescope" with a resolution of 
~5 x 10~ 7 arcseconds (given an angular diameter distance to the quasar of ~ 5.4 x 10 27 cm). 
In other words, in SDSS J1004+4112, nature has provided us with an extremely powerful 
tool for the study of BELR structure. 
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3.1.4- Some Predictions 

The microlensing hypothesis leads to several predictions that may guide further obser- 
vations. First, the nature of the object(s) responsible for the microlensing is unknown. Two 
obvious possibilities are stars or MACHOs associated with a galaxy (presumably a cluster 
member) with a small projected impact parameters to component A, or stars or MACHOs 
in the intracluster medium. While the Subaru i-band image of the field (Fig. 3) suggests 
that there is indeed a galaxy near component A, the saturation of the quasar image makes 
PSF subtraction uncertain. Deep, high-resolution images, preferably in the near-IR, would 
better reveal whether this galaxy is real. If so, then the microlensing optical depth would 
be relatively high for component A and (perhaps quite) low for the other components, so 
further monitoring would probably reveal additional microlensing events in component A 
but not in the other components. 

Alternatively, if the microlensing is caused by intracluster stars or MACHOs then the 
microlensing optical depth for component A is probably fairly low — but it is likely to be 
similar for the other components. In this case further monitoring could well reveal microlens- 
ing events in the other components, and the frequency of events would reveal the number 
density of intracluster objects. 

If additional microlensing events are detected, it will be extremely interesting to track 
which lines are microlensed as a function of time. That would provide a unique and powerful 
new probe of ionization stratification in the BELR. In addition, the velocity dependence of 
any enhancements of the BELR in future microlensing events will help to reveal the kinematic 
structure of the BELR. 



3.2. Time-independent Emission Line Differences 

In addition to the time-dependent differences in the blue wings of the high-ionization 
emission lines, there are also subtle differences in the C IV emission line profiles between the 
components that persist over at least 322 days in the observed frame. In Figure 2 we have 
overplotted a Gaussian at the position of C IV emission to guide the eye and help illustrate 
these differences. 

In particular, we note that there is a slight excess of high-velocity redshifted C IV 
emission in component B as compared with component A. This excess is best seen in the 
2003 May 31 Keck spectra, where there is a kink in the profile of B near A 1560 A that is 
not present in component A. In addition, in all observations of component B the fall-off in 
the red wing toward A 1600 A is more gradual than in component A. We also see similar 
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(but opposite) differences in the blue wing with the most extreme blue wing flux falling off 
more gradually in component A than in component B. In other words, at all epochs the 
C IV emission line profile is somewhat blueward asymmetric in component A, and somewhat 
redward asymmetric in component B. 

The cause of these differences is unclear. Because the predicted A-B time delay is <30 
days, they are unlikely to be due to intrinsic variability (unless the time delays are grossly 
underestimated). One interesting possibility is that they are also due to microlensing, but 
with a mass scale that is much larger than for the blue-wing BELR microlensing discussed 
above (in order to make the variability time scale longer than 322 days in the observed 
frame). The responsible objects could perhaps be globular clusters, or clumps of dark matter 
of mass ~10 4 -10 8 M Q , in which case the phenomenon would be termed millilensing (rather 
than microlensing) and could provide a unique probe of dark matter substructure of the type 
predicted in the Cold Dark Matter paradigm (Metcalf & Madau 2001; Dalai & Kochanek 
2002; Wisotzki et al. 2003; Metcalf et al. 2003). 

One way to test the millilensing hypothesis would be to normalize the spectra using 
narrow lines since they should be insensitive to small-scale structure (Moustakas & Metcalf 
2003). Any differences in the broad lines would then indicate millilensing. This might be 
possible with either J-band IR spectra of [O III] AA4959, 5007 or with the unusually strong 
(but not apparently microlensed) nitrogen lines seen in our optical spectra, since the nitrogen 
lines appear relatively narrow and may lack a broad component. 

4. Conclusions 

We have presented seven epochs of spectroscopic data on the two brightest components 
of the wide-separation, quadruply imaged quasar SDSS J1004+4112. Although the simplest 
lensing scenarios predict that the four components should have identical spectra, the data 
reveal significant differences in the emission line profiles of the components. In particular, 
the C IV emission line profile in components A and B show both variable differences and 
differences that are constant over 322 observed-frame days. The He II and Si IV/O IV] lines 
in component A also show variability similar to that seen in the C IV line. 

Because the predicted time delay between A and B is <30 days, we argue that the 
differences are not due to intrinsic variability in the quasar coupled with a lensing time delay. 
Instead, we suggest that the variability in the blue wing of component A is best explained by 
microlensing of part of the broad emission line region, resolving BELR structure on the order 
of a few light days. This represents the first robust detection of BELR microlensing, with 
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evidence based on multiple emission lines and involving observed variability. The nature of 
the time-independent differences is less clear, but they may also be the result of a lensing 
event. In any case, it is clear that continued spectroscopic monitoring of SDSS J1004+4112 
should be carried out in an attempt to map the structure of its broad emission line region 
through additional microlensing events. 

Funding for the creation and distribution of the SDSS Archive has been provided by 
the Alfred P. Sloan Foundation, the Participating Institutions, the National Aeronautics 
and Space Administration, the National Science Foundation, the U.S. Department of En- 
ergy, the Japanese Monbukagakusho, and the Max Planck Society. The SDSS Web site is 
http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium 
(ARC) for the Participating Institutions. The Participating Institutions are The University 
of Chicago, Fermilab, the Institute for Advanced Study, the Japan Participation Group, The 
Johns Hopkins University, Los Alamos National Laboratory, the Max-Planck-Institute for 
Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State 
University, University of Pittsburgh, Princeton University, the United States Naval Obser- 
vatory, and the University of Washington. 

Based on observations obtained with the Apache Point Observatory 3.5-meter telescope, 
which is owned and operated by the Astrophysical Research Consortium. Based on observa- 
tions obtained at the W. M. Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of California and the National 
Aeronautics and Space Administration, made possible by the generous financial support of 
the W. M. Keck Foundation. Based in part on data collected at Subaru Telescope, which is 
operated by the National Astronomical Observatory of Japan. G. T. R. was supported in 
part by HST grant HST-GO-09472.01-A. C. R. K. is supported by NASA through Hubble 
Fellowship grant HST-HF-01141.01-A from the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. 



-13- 



REFERENCES 

Abajas, C, Mediavilla, E., Muhoz, J. A., Popovic, L. G, & Oscoz, A. 2002, ApJ, 576, 640 

Arav, N., Barlow, T. A., Laor, A., Sargent, W. L. W., & Blandford, R. D. 1998, MNRAS, 
297, 990 

Baltz et al. 2003, astroph/03 10845 

Cao, L., Wei, J.-Y., & Hu, J.-Y. 1999, A&AS, 135, 243 

Chartas et al. 2003, astroph/0401240 

Colley, W. N. & Schild, R. E. 2003, ApJ, 594, 97 

Dalai, N. & Kochanek, C. S. 2002, ApJ, 572, 25 

Elvis, M. 2000, ApJ, 545, 63 

Eracleous, M. 2003, astroph/0310649 

Eracleous, M. & Halpern, J. R 2003, ApJ, 599, 886 

Filippenko, A. V. 1989, ApJ, 338, L49 

Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shimasaku, K., & Schneider, D. P. 1996, 
AJ, 111, 1748 

Gunn, J. E., Carr, M., Rockosi, C, Sekiguchi, M., Berry, K., Elms, B., de Haas, E., Ivezic , 
Z., et al. 1998, AJ, 116, 3040 

Huchra, J., Gorenstein, M., Kent, S., Shapiro, I., Smith, G., Horine, E., & Perley, R. 1985, 
AJ, 90, 691 

Inada et al. 2003, Nature, 426, 810 

Irwin, M. J., Webster, R. L., Hewett, P. C, Corrigan, R. T., & Jedrzejewski, R. I. 1989, AJ, 
98, 1989 

Kaspi, S., Smith, P. S., Netzer, H., Maoz, D., Jannuzi, B. T., & Giveon, U. 2000, ApJ, 533, 
631 

Kayser, R., Refsdal, S., k Stabell, R. 1986, A&A, 166, 36 
Kochanek, C. S. 2003, astroph/0307422 



-14- 



Lewis, G. F. & Ibata, R. A. 2003, MNRAS, astroph/0310818 
Metcalf, R. B. & Madau, P. 2001, ApJ, 563, 9 
Metcalf et al. 2003, astroph/0309738 

Moustakas, L. A. & Metcalf, R. B. 2003, MNRAS, 339, 607 
Murray, N. & Chiang, J. 1998, ApJ, 494, 125 
Nemiroff, R. J. 1988, ApJ, 335, 593 
Oguri et al. 2003, ApJ, astroph/0312429 

Oke, J. B., Cohen, J. C, Carr, M., Cromer, J., Dingizian, A., Harris, F. H., Labrecque, S., 
Lucinio, R., et al. 1995, PASP, 107, 375 

Peterson, B. M. & Wandel, A. 1999, ApJ, 521, L95 

Richards, G. T., Fan, X., Newberg, H. J., Strauss, M. A., Vanden Berk, D. E., Schneider, 
D. P., Yanny, B., Boucher, A., et al. 2002, AJ, 123, 2945 

Schechter, P. L., Udalski, A., Szymahski, M., Kubiak, M., Pietrzyhski, G., Soszyhski, I., 
Wozniak, P., Zebruh, K., et al. 2003, ApJ, 584, 657 

Schechter, P. L. & Wambsganss, J. 2002, ApJ, 580, 685 

Schneider, P., Ehlers, J., & Falco, E. E. 1992, Gravitational Lenses ( Springer- Ver lag: Berlin) 
Schneider, P. & Wambsganss, J. 1990, A&A, 237, 42 

Strateva, I. V., Strauss, M. A., Hao, L., Schlegel, D. J., Hall, P. B., Gunn, J. E., Li, L., 
Ivezic, Z., et al. 2003, AJ, 126, 1720 

Totani, T. 2003, ApJ, 586, 735 

Wambsganss, J., Paczynski, B., & Katz, N. 1990, ApJ, 352, 407 
Wandel, A., Peterson, B. M., & Malkan, M. A. 1999, ApJ, 526, 579 

Wisotzki, L., Becker, T., Christensen, L., Helms, A., Jahnke, K., Kelz, A., Roth, M. M., & 
Sanchez, S. F. 2003, A&A, 408, 455 

Wozniak, P. R., Udalski, A., Szymahski, M., Kubiak, M., Pietrzyhski, G., Soszyhski, I., & 
Zebruh, K. 2000, ApJ, 540, L65 



- 15 - 

York, D. G., Adelman, J., Anderson, J. E., Anderson, S. F., Annis, J., Bahcall, N. A., 
Bakken, J. A., Barkhouser, R., et al. 2000, AJ, 120, 1579 



This preprint was prepared with the AAS IATgX macros v5.2. 



-16- 
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Wavelength (Angstroms) 

Fig. 1.— Keck/LRIS spectra of components A, B, C, and D of SDSS J1004+4112. Dashed 
vertical lines indicate the expected peaks of Si IV/O IV], C IV, and He II for z = 1.734. A 
power law continuum, fit between A 1450 A and A 1690 A, has been subtracted from each 
spectrum, and the spectra are all normalized to the peak of C IV. The spectra are smoothed 
by a seven pixel boxcar filter. 
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Fig. 2. — Seven epochs of data in the C IV emission line region of SDSS J1004+4112. A power 
law continuum, fit between A 1450 A and A 1690 A, has been subtracted from each spectrum, 
and the spectra are all normalized to the peak of C IV. We also overplot a scaled Gaussian at 
the center of C IV to guide the eye toward emission line differences that are persistent with 
time. The dashed vertical line indicates the expected C IV peak at A 1549.06 A. All spectra 
have been smoothed to similar resolutions. Note that this smoothing hides the associated 
absorption system that is observed just blueward of the C IV emission line peak. 
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Fig. 3. — Central (22'.'6) region of the 1340 s i-band Subaru Prime Focus Camera image from 
Oguri et al. (2003, Fig. 8). North is up, East is left. Although the image is saturated and 
proper PSF subtraction of component A is not possible, there appears to be a superimposed 
galaxy (with an estimated % ~ 24.5) just to the North-East of component A that could 
host the microlensing object. The white circle is centered on the possible microlensing 
galaxy and represents a 20/if 1 kpc radius at z = 0.68 assuming VLm = 0.3, Vt^ = 0.7, and 
H = TO/irokms^Mpc- 1 . 
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Fig. 4. — Comparison of Keck spectra of components A and B in the C III] to Mg II region 
of the spectra, showing the weaker enhancement (relative to the higher ionization lines) of 
the blue wings of these emission lines. A power law continuum, fit between A 1690 A and 
A 2200 A, has been subtracted from each spectrum, and the spectra are all normalized to 
the peak of C III]. No constraint is placed on Mg II. The dashed vertical lines indicate the 
expected Al III, C III] and Mg II peaks in addition to the Fe III UV34 and UV48 complexes. 
The spectra are smoothed by a seven pixel boxcar filter. 
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Table 1. Summary of Observations 



Date (UT) 


Telescope 


Components 


Exp. Time 


Dispersion 


A Range 








(sec) 


(Apix- 1 ) 


(A) 


2003 February 3 


SDSS 2.5m 


B 


2700 


~ 1 


3800-9200 


2003 May 3 


ARC 3.5m (DISIII) 


A 


1800 


~ 2.4 


3820-5630 


2003 May 31 


Keck I (LRIS) 


ABCD 


900 


~ 1 


3028-9700 


2003 November 21 


ARC 3.5m (DISIII) 


AB 


2700+2700 


2.4 


3890-9350 


2003 November 30 


ARC 3.5m (DISIII) 


AB 


2700+2700 


2.4 


3890-9350 


2003 December 1 


ARC 3.5m (DISIII) 


AB 


2400+2400 


2.4 


3890-9350 


2003 December 22 


ARC 3.5m (DISIII) 


AB 


2400+2400 


2.4 


3890-9350 



